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ABSTRACT 


The behavior of spatially growing three-dimensional waves in a supersonic boundary 
layer is studied. The objectives are to validate our computer code by comparing with ex- 
periments and parallel and non-parallel linear stability theories for = 4.5, yjRe/ft = 
1550 and at T ti00 = 100°F and then use the program to study nonlinear effects. The 
three-dimensional unsteady Navier- Stokes equations are solved by a finite difference 
method which is fourth-order and second-order accurate in the convection and viscous 
terms respectively, and second-order accurate in time. Spanwise periodicity is assumed. 
The inflow disturbance is composed of eigenfunctions from linear stability theory. Com- 
puted results for small amplitudes of this disturbance agree well with linear theory and 
experiment for several frequencies. By increasing the amplitude of the inflow distur- 
bance, nonlinear effects in the form of a relaxation type oscillation of the time signal of 
pu are observed. 


1 This research was supported in part by the National Aeronautics and Space Administration under 
NASA Contract Nos. NAS1-18107 and NAS1-18605 while the second author was in residence at the 
Institute for Computer Applications in Science and Engineering (ICASE), NASA Langley Research 
Center, Hampton, VA 23665. 
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1. INTRODUCTION 


The objective of this paper is to study the linear and nonlinear behavior of spatially 
growing waves in a supersonic flat plate boundary layer. The behavior of these waves 
will be obtained by numerically solving the three-dimensional, unsteady Navier-Stokes 
equations using a fourth-order accurate finite difference scheme. The numerical results 
will demonstrate a linear regime where the growth is well predicted by linear stabil- 
ity theory and a nonlinear regime characterized by the generation of higher harmonics 
associated with the appearance of relaxation type oscillations. 

The numerical computations will be conducted using a two-dimensional mean flow 
over a flat plate. The parameters of the mean flow will be chosen to agree with ex- 
perimental measurements (Kendall [1967]). The mean flow will be perturbed by a dis- 
turbance at the inflow boundary. This disturbance will be taken as an eigenfunction 
obtained from linear stability theory (El Hady [1981]). Due to the spatial growth of the 
waves, nonlinear behavior can be obtained in two ways: either by increasing the ampli- 
tude of the disturbance at inflow or by following the disturbance for a longer distance 
downstream. Both mechanisms will be noted here. 

An extensive discussion of linear stability theory can be found in Mack [1969]. The 
results in Mack [1969] were obtained by assuming parallel flow and solving the Orr- 
Sommerfeld equations for the vertically varying eigenfunction. Nonparallel effects can 
be modeled by allowing the growth rate to vary with the longitudinal coordinate x. 
Nonlinear effects cannot be accounted for in this theory. The growth rates predicted 
by linear stability theory were compared with experimentally measured growth rates of 
Kendall [1967] and found to be in fair agreement for a flow at Mach number 4.5 (the 
mean flow used in the numerical computations). An approximation to the nonparallel 
effect of the same Mach number has been investigated by El Hady [1981] using a multiple 
scale expansion. His results are also in fair agreement with Kendall’s experiment and 
Mack’s analysis. 

Nonlinear spatial instability was studied for subsonic flows in two dimensions by 
Bayliss, et al. [1985, 1968], in three dimensions by Parikh et al [1987], and for incom- 
pressible flows in three dimensions by Fasel, et al. [1987] and Murdock [1986]. The 
computations in Bayliss, et al. [1986] demonstrated a progressive sharpening of the time 
signals as the disturbance propagates downstream. At a sufficiently large downstream 
location, the fluctuating disturbances take on the character of relaxation oscillations. 
Similar results are shown here for the three- dimensional computations in supersonic 
flow. The computations also show a progressive clustering of vorticity contours as the 
wave propagates downstream and an uplifting of vorticity away from the wall. Ultimately 
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the nonlinear growth will saturate and the flow will undergo transition. 

In this paper we study the nonlinear distortion, occurring prior to transition, with 
increasing spatial distance. Computations of temporally growing waves assuming pe- 
riodicity in the streamwise direction have been reported in Erlebacher and Hussaini 
[1986]. The study of spatially growing instability waves as conducted here will permit 
the study of nonparallel effects on the nonlinear growth. In addition, the growth of the 
disturbance will be generated in a more physically realizable manner. The additional 
resolution requirements needed to resolve the spatial growth are severe and this limits 
the degree of nonlinearity that can be reliably computed. In Section 2, we describe the 
numerical procedure, and in Section 3, we discuss our numerical results. 


2. COMPUTATIONAL METHOD 

The numerical method uses an explicit finite difference scheme to solve the unsteady 
Navier-Stokes equations in conservation form. The finite difference scheme is fourth 
order accurate in the spatial coordinates and second order accurate in time. The scheme 
is a generalization of MacCormack’s scheme due to Gottlieb and Turkel [1976] and has 
been previously used in the study of two-dimensional spatially unstable waves (Bayliss, 
et al. [1985, 1986]). 

The Navier-Stokes equations can be written as 

w t = F x + G y + H z (2.1) 

where w is the vector (p, pu, pv , pw , E) T , p is the density, u, v, and w are the streamwise, 
vertical, and spanwise components of the velocity vector and E is the internal energy. 
The functional form of the flux functions F, G, and H are well-known and for the sake 
of brevity will not be given here. For the scalar equation 


u t = fy 


we consider the predictor-corrector scheme to compute 


u 


,n+l 


~ u((n + l)At, jx) 


Uj - U" + + 8/j+i - /j+2 ] 


( 2 . 2 ) 


u i +1 = 2 + A* + g2y^[Ui ” + U- 2 ]] 

where fj = It is known that if (2.2) is alternated with its symmetric variant, the 

resulting scheme is fourth order accurate provided At = 0(Ax 2 ). For the case of flux 
functions which depend on derivatives of u, (e.g., as in (2.2)), the scheme is fourth order 
accurate on the inviscid terms and second order accurate on the viscous terms. 
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In the numerical method used here, the scheme (2.2) is implemented with opera- 
tor splitting, Strang [1968]. That is, if L x denoted the solution (via 2.2)) to the one- 
dimensional equation 

wt + fx = 0 (2.3) 


and L y ,L z are defined similarly, then the solution is advanced two timesteps via the 
formula 


w n+2 = BL x L v L z BL z L v L x w n 


(2.4) 


where the operator B denotes the boundary conditions (discussed below). The resulting 
scheme is fourth order accurate provided the timestep is sufficiently small. A compar- 
ison of this scheme with the second order MacCormack schemes for two-dimensional, 
spatially-unstable waves was presented in Bayliss, et al. [1985] and indicated a signifi- 
cant improvement in accuracy for the higher order scheme. 

The solution is obtained in a rectangular box illustrated by Figure 1. Spanwise 
periodicity is assumed. At the wall, the no-slip condition and a wall temperature T w are 
imposed, i.e., 

u = v = 0 


T = T W . 


(2.5) 


At the inflow boundary, the field is imposed. The inflow data is of the form 


w = wq + s Real Part [p*^ t+ ^*^(y)]. (2.6) 


Here, w 0 is the mean state, a spreading two-dimensional boundary layer computed in 
a previous calculation. The eigenfunction (( y ) is obtained from the stability code of 
El Hady [1981] and the frequency u> and spanwise wave number (3 are specified. The 
parameter e measures the strength of the inflow disturbance. At the outflow boundary, 
the solution is extrapolated from the interior. At both boundaries there is a small 
subsonic region near the wall where one less boundary condition should be imposed at the 
inflow and one boundary condition should be imposed at the outflow. These segments are 
treated by using one-dimensional normal characteristics described in Bayliss and Turkel 
[1982] and similar to the treatment of the upper boundary described immediately below. 

The upper boundary is treated as a subsonic outflow boundary as v is small and 
positive. At this boundary, the outgoing normal characteristic variable ( p — pcv) based 
on linearizing the inviscid equations about the current state is prescribed as its initial 
value. The remaining characteristic variables are extrapolated from the interior. This 
approach was used in Bayliss and Turkel et al. [1982] and found to be effective in 
allowing the upper boundary to be brought closer to the edge of the boundary layer. 
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For the computations presented here, the vertical extent of the computational domain 
is approximately 5 Si where Si is the boundary layer thickness at the inflow. 

All of the computations presented here were obtained with a mean flow with freestream 
Mach number M m = 4.5, total temperature T too = 100°F and Reynolds number per foot, 
Re/ft = 2.425 x 10 6 . This corresponded to the experimental conditions of Kendall. The 
inflow boundary was one foot from the leading edge. The Reynolds number based on 
boundary layer thickness Re$* was 17900. The wall temperature T w was 110. 8°R. 

3. NUMERICAL RESULTS AND DISCUSSION 

We first consider inflow disturbances of small amplitude so as to remain in the linear 
regime. This enables a study of complete nonparallel effects and a comparison or our 
results to linear stability theory and also to the experimental measurements of Kendall 
[1967]. 

We first consider an inflow disturbance level of .075% of the free stream. One measure 
of the solution is the growth of the disturbance downstream. We compute our growth 
from the fluctuating mass flux 

(pu)' = pu- < pu> (3.1) 

where < pu > is the computed mean in time. 

In order to quantify the growth of the disturbance, we compute the root mean square 
(RMS) in t of (pu)'. This is then integrated in y and z, i.e., 

G(x ) = J (RMS(pu)' )dzdy . (3.2) 

The growth rate o is then obtained from an approximation to l/xlog(Cr(x)/G(xo)). 
We have computed cr(F) for a range of non-dimensional frequencies F. Two modes, the 
first and second modes are of interest. The first mode is three-dimensional and for the 
parameters considered here the most amplified wave propagates at an angle of 55°. It is 
known, Mack [1969], that for the parameters considered here the most unstable second 
mode is two dimensional. The most amplified second mode has an amplification rate 
approximately three times greater than the most amplified first mode. 

In Figure 2, we indicate the computed values of growth rate c(F) for several fre- 
quencies and for both first and second modes. The data is compared to the growth rates 
predicted from linear theory (Mack [1969]), experiment (Kendall [1967]) and nonparallel 
theory (El Hady [1981]). In all cases, the peak inflow disturbance was .075% of the free 
stream (in terms of u). We note that both the nonparallel results and the experimental 


5 


data are for the growth of the maximum disturbance, not for the integrated R.M.S. 
disturbance as in (3.2). 

The numerical results show that for all frequencies considered the computed growth 
rate lies close to both the parallel and nonparallel results. Where experimental results 
are available, the computed growth rates are close to the experimental measurements. 
The results show that the magnitude of full nonparallel effects (i.e., as obtained from 
the Navier-Stokes equation) are of the same order as that predicted from nonparallel 
theory. 

The computational domain extends approximately 8.3 wavelengths in the streamwise 
direction (based on the value of A x at the inflow determined from linear theory. There 
are 28 points per wavelength in the streamwise direction and 150 points in the vertical 
direction (using an exponential stretching to concentrate points in the boundary layer). 
When £ = 0.075 % and 2%, 15 points per wavelength in the spanwise direction were 
used. When e = 3%, 32 points per wavelength were used. A case with e = 1.5% 
was run with 15 and 30 points per wavelength in the spanwise direction and virtually 
identical results were obtained in the early stages, i.e., before pronounced nonlinear 
effects became apparent. However, the spatial resolution in the spanwise direction, even 
with 32 points, will not be sufficient when nonlinear effects become pronounced as in 
the case with e = 3%. An effort to change the spanwise approximation to one using a 
Fourier spectral method is now underway. 

In Figure 3, we present the time signal of pu(t) (non-dimensional) at a fixed location 
in the boundary layer and for several downstream stations. This is for the case where 
the inflow disturbance is 3% and nonlinear effects should be expected . Close to the 
inflow boundary, the oscillations are sinusoidal. With increasing distance downstream, 
weak nonlinear effects first emerge in the form of relaxation oscillations. There is a 
broadening of the valleys of the wave form and a narrowing of the peaks. (The opposite 
pattern was noted at other points in the boundary layer). At x = 5A X from the inflow 
(Re^ = 2.3xl0 4 ) what may be the initial stages of breakdown are noticeable. At a lower 
value of the inflow perturbation, these effects are not as pronounced. In Figure 4, for 
e = 1%, at 3A X from the inflow, pu(t) at the same location in the boundary layer as 
in Figure 3 is plotted. The oscillations are not of the relaxation type. In Figure 5, 
spanwise- vorticity contours are plotted. In these figures Si denotes the inflow boundary 
layer thickness. This is for the case with e = 3% and F = 1.5xl0 -5 and at a time 
t = 0.89. Note that the disturbance has not reached the outflow boundary. Figure 6 
has two such plots at an earlier time. They do not show as advanced a state. Figure 
7, where t = 0.99 is close to the point where numerical instability prevented us from 
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continuing our computations any further. It shows a much more advanced state than 
t = 0.89. The spatial resolution here is suspect. 

Several physical features of a boundary layer approaching transition are noticed in 
Figure 5. There is a lifting of vorticity away from the wall towards the edge of the bound- 
ary layer. Also, the bending downward of these vortices is evident. This phenomenon, 
the one spike stage, is one of the early stages of transition. 

In Figure 8, contours of vorticity in the x — z plane are presented. Streamwise 
vorticity contours at three downstream locations are shown in Figure 8. Both these 
sets of figures also show the progresion (with increasing downstream distance) towards 
breakdown. The time development of this process is also apparent from Figure 9. 

Finally in Figure 10, the instantaneous velocity on the plane z = are shown. The 
perturbation imposed at the inflow only modifies the profile slightly but the distortion 
becomes more pronounced with increasing distance downstream. At a downstream lo- 
cation where nonlinear effects are apparent from spanwise vorticity plots, i.e., between 
4A S and 5A X (Re$ between 2.2 x 10 4 and 2.3 x 10 4 ), the velocity profiles indicate a strong 
distortion. 

4. CONCLUSION 

The three-dimensional stability of a growing, supersonic flat-plate boundary layer 
was studied numerically by solving the complete Navier-Stokes equations. Satisfactory 
comparison with linear parallel and non-parallel stability theories, and experiment are 
obtained when a small amplitude inflow disturbance is used. For higher inflow dis- 
turbance amplitudes, nonlinear effects are apparent. The nonlinear effects obtained in 
the present study are only the early stages of transition. Several modifications to im- 
prove the accuracy of the code will be necessary before it can reliably compute the 
more advanced stages of transition. These would include a variable grid spacing in the 
streamwise direction and a spectral approximation in the spanwise direction. 

Work in progress includes the incorporation of some of these changes as well as the 
study of the nonlinear interaction of two- and three-dimensional waves. 
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Schematic of three-dimensional computa- 
tional domain used. 
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Figure 2. Comparison of computed growth rates 

with theory and experiment. 
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Figure 3. 
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Time series of p u at | = 0.795 for e = 3% 
case. This shows the development of re- 
laxation oscillations with distance down r 
stream of inflow. 
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Figure 4. Time series of pu at | = 0.795 for e = 

1%. Comparison with fig. 4 shows that re- 
laxation oscillations are not as pronounced 
in this case. 
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Development of vorticity contours in the x — z 
plane at y/6i = 0.795. 
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Instantaneous velocity profiles on the f- = 
| plane. 
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